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>t ABSTRASPMOANOND 


An experimental stainless steel heat pipe using water as the working 
fluid and 400 mesh stainless steel screen for a wick was designed and 
tested to determine the effect of gravity and nucleate boiling on heat 
pipe performance. The results of heat pipe operation at various angles 
of inclination in a gravity field are presented and compared with the 
existing theoretical predictions. 

The maximum heat flux obtained experimentally at angles of inclina- 
tion less than 90 degrees was less than the predicted value by a factor 
of two or three. The maximum heat flux obtained for an angle of inclina- 
tion greater than 90 degrees was much higher than that predicted. In 
addition, nucleate boiling noise was detected during operation at 
angles of inclination greater than 90 degrees. 

The experimental results coupled with visual examination of the 
pipe after operation indicate that the pipe was not performing satis- 
factorily. Recommendations for a better design of an experimental heat 


pipe are presented. 
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NOMENCLATURE 
Area, fa 
Screen permeability factor 
Acceleration due to gravity, et Mec” 
Latent heat of vaporization, BTU/1b 
Mass rate of flow, lb/sec 
Heat flow, watts, BTU/sec, BTU/hr 
Radius of vapor passage, ft 
Reynolds number = Q/7T LRN 
Outside radius of wick structure, ft 
Effective capillary radius of wick, ft 
Heat pipe length, ft 
Surface tension, lb/ft 
Incremental error 
Wick void fraction 
Dynamic viscosity of the vapor 
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Kinematic viscosity 
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I. INTRODUCTION 

Background 

The transportation of large amounts of heat plays an increasingly 
important role in present day engineering problems and applications. 
It is not surprising that a major effort has been devoted to research to 
find more efficient ways to transfer large amounts of heat. The heat 
pipe is one such device that has recently been developed which is capable 
of transferring large amounts of heat efficiently. Copper, which is 
usually regarded as one of the best conductors, is not nearly as efficient 
as the heat pipe. For instance, if thermal power of 10,000 watts were 
applied to one end of a solid copper bar one inch in diameter and one foot 
long, the temperature difference along the bar could theoretically exceed 
30,000°F. This same amount of thermal power could be transferred through 
a heat pipe of similar dimensions with only a small temperature difference 
[1]. In this Heine wabene is transferred by the evaporation of a liquid 
in one section of the pipe and the condensation of the vapor at a colder 
section of the pipe. An ordinary reflux capsule (or thermosiphon) trans- 
fers heat in this manner, also. The unique property of the heat pipe is 
that no outside source of energy is required to return the condensate to 
the evaporator. In a heat pipe, the return of the condensate is accomplished 
by the use of porous wicking material (such as wire screen, small channels, 
or sintered porous material) inside the pipe. The heat pipe is a self- 
contained, closed system requiring no outside source of energy other than 
the energy being transported. 

The basic principle of the heat pipe was first discovered by R. S. 
Gaugler in 1942 but was not utilized until 1963, when G. M. Grover inde- 


pendently discovered the device and named it a "heat pipe" [1]. 


Since 1963, Grover [2], Cotter [3,4], and others have done research, both 
theoretically and experimentally, to determine design and operating 


characteristics of these devices. 


Theory of Operation 


The heat pipe is essentially a closed, hollow structure whose inside 
wall is covered by a porous wick. It is not necessary that the wick be 
placed against the inside surface of the container, although this will 
normally be the best place for it. Condensation and evaporation take 
place at the lLiquid-vapor interface, therefore if the wick is placed 
against the inside surface the radial transfer of heat necessary for 
operation will occur through the medium of highest thermal conductance 
and radial temperature differences will be minimized. If for some reason 
it is not desirable to place the wick next to the wall, it is possible to 
use fins to conduct the heat to the wick and thus reduce the radial tempera- 
ture difference. 

For theoretical simplicity, assume the heat pipe is a long, hollow, 
closed cylinder with the wick adjacent to the inside surface. (Figure 1). 
The theory for steady state operation for such a pipe has been developed 
by Cotter [3] and is briefly stated below. 

The assumptions made in Cotter's theory are: 

(1) The pipe has a large length to diameter ratio. 

(2) The pipe contains sufficient pure liquid to fill the wick. 

(3) The vapor space contains no non-condensable gases. 

(4) The vapor flow is in the incompressible regime and inertial 
forces predominate. 

(5) The pressure distribution in the wick obeys the hydrostatic 
laws for an incompressible fluid. 


(6) Flow through the wick may be given by a form of Darcy's law 
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for flow through porous media. 

(7) The diameter of the pipe is constant. 

(8) Heat addition and removal is uniform. 

(9) Nucleate boiling does not occur in the wick. 
Axial heat transfer is accomplished by flow of the vapor from the hot 
end of the pipe to the cold end. The vapor is condensed at the cold 
end of the pipe and is returned to the hot end by means of the wicking 
material. The driving force for successful operation of the pipe is 
capillary action. 

In order to arrive at the equations governing the operation of the 


heat pipe, consider the sketch below. The sketch represents a partial 





cross-section of the heat pipe which schematically shows the wick and 
part of the vapor cavity. The vapor pressure at a given axial position 
is given by P tz). Likewise, the liquid pressure at some axial position 
is given by Piz). The subscripts c and e represent positions in the 
condenser and evaporator sections, respectively. Writing a pressure 


balance for the pipe gives the following expression: 
= 2s E be 
(P(2,)-P (2+ GRY Pa + GB (2 Der @.)) 
+t. a = 
(P, (2,)-P (z,)) = 0 
Now at any axial position, z, the interface between the liquid in the wick 


and the adjacent vapor must assume a local radius of curvature, r(z), such 


TZ 


that the surface tenxion, if , supports the difference in pressure between 


the liquid and the vapor. Thus, 


B(2) Pa aaae § cos 6/r(z) 
where 9 is the contact angle between the surface of the liquid and the 
capillary structure. In the condenser, the radius of curvature, r(z_); 
will be large due to condensation and may approach infinity while the 


contact angle 9 approaches 90 degrees so that 
= = Ne 
Patz.) a2, (ay). 22 S cos @/r(z_) 0 
On the contrary,in the evaporator section, the radius of curvature, r(z.); 


will be small due to evaporation and 9 approaches zero so that the pres- 


sure drop across the liquid-vapor interface may be written 
P (zd -P, (2) = 2 o /r(z,) 
Letting P(z,)-P(z_) = Aap and P, (2,)-P, (z,) =Br,, 


with the above simplifications, then the pressure balance may now be re- 


written as 
L\ P + AP, = 2 F /r(z,) (1) 


Darcy's law for flow through porous media yields an expression for 


the pressure drop in the liquid 


Ap “ by ZQ ‘ ‘1 gE cop. ® (2) 
1 4 Z g 
Z IT (Ro RY Er oL O 


where b is a dimensionless constant which depends on the geometry of the 
porous wick, \J is the kinematic viscosity, Z is the length of the pipe, 

Q is the heat applied, R is the radius of the vapor cavity, and R is the 
outer radius of the wick. The wick void fraction is represented by € and 
rn. is the effective capillary radius. L is the latent heat of vaporization 


of the fluid, Py is the density of the liquid, g/g, is the gravitation 
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term and @ is the angle of inclination of the pipe with the vertical. 

The effective capillary radius, re for a wire screen or porous 
wicking material must be determined by a capillary rise experiment in 
which the height of liquid supported by the wicking material is experi- 
mentally measured. This height h, is used in the relationship 

R= 2¢c08 6/69/90) h 
to define ene Wereative capillary radius. In the above relationships, /(O 
is the density of the liquid used in the experiment, and cos @ is assumed 
to be unity. If the capillary structure is a small tube of circular cross- 
section, then the effective capillary radius turns out to be the inner 
radius of the tube as would be expected. 

Two major regimes exist for flow in the vapor phase. These regimes 
depend on the magnitude of a Reynolds number, Ro based on the radial flow 
velocity of the vapor at the channel wall, the vapor cavity radius, R, the 
vapor density, Py and dynamic viscosity, i The previously assumed 
conditions for vapor flow yield the following relationship for the pres- 
sure drop in the vapor phase: 


4, \ 22 : R. Z< 1 


A Bos 3 2 2 (3) 


The above equations represent the steady state solutions for the pressure 
drops in the liquid and vapor phases. In general, the Reynolds number, R.» 
will be greater than l. 

The local radius of curvature, r(z), will be at least as large as the 
vapor cavity radius, R, and may be as small as the effective capillary 
radius of the wick, ro: The effective capillary radius determines the 


maximum pressure difference which may exist between the liquid and the 
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vapor. If r(z) is the effective capillary radius, te. equation (1) may 


be rewritten as 
BGs +) flee, geo fe (4) 


Substituting equations (2) and (3) into (4) yields 


(Tete OF Ayes b V 2Q +e coe Bn 
8 p RL 2 TY (R*-R*) E rb B, ey: 
(5) 


The solution of equation (5) for Q gives the maximum heat flux that can be 
attained in a heat pipe which is limited by capillary action only. 

Physical indications of having reached the maximum heat transfer limit 
are a sudden rise in temperature in the evaporator and a sudden drop in 
temperature in the condenser. 

It can be seen from equation (5) that a number of factors govern the 
operation of a heat pipe. These factors may be grouped into three broad 
categories: the maximum heat flux, the pipe geometry and size, and the 
properties of the working fluid. By specifying the proposed application 
and selecting a working fluid for the temperature range desired, one can 
design a heat pipe of size and geometry to fit the application. In other 
words, the application of the pipe and the working fluid determine the 
design of the pipe. 

In deriving the equations for the operation of a heat pipe, the as- 
sumption was made that nucleate boiling did not occur. The reason for the 
above assumption was because it was presumed that the formation of bubbles 
in the wick would prevent the condensate from returning to the evaporator 
or would at least impede flow to the evaporator. The result of insuffi- 
cient flow to the evaporator would be the formation of hot spots and pre- 


mature drying out of the evaporator. 


1 


Previous Investigations 


Since the work by Grover in 1963, much work has been done investi- 
gating the properties and uses of the heat pipe. In 1965, Grover, et. al. 
[2] investigated the use of a heat pipe as a heat removal device for 
thermionic devices in space. B. D. Marcus [5], apparently working independ- 
ently at TRW Space Technology, arrived at a theory to determine the maximum 
heat flux for a heat pipe. Marcus arrived at a method which optimized the 
thickness of the wick such that the wick was just large enough to permit 
the condensate flow rate corresponding to the maximum heat transfer rate 
and just small enough to avoid nucleate boiling in the wick at the maxi- 
mum heat transfer rate. Recently Cotter [4] investigated the transient 
modes of heat pipe operation. He discussed flow properties during startup 
and the mechanisms by which startup can fail. Two phenomena can cause 
startup to fail. The first one being the attainment of the superheating 
limit at which hot spots form in the evaporator excluding the liquid from 
where it is needed due to boiling in the capillary structure. The second 
includes all other mechanisms of failure which result in insufficient re- 
turn of condensate to the evaporator. J. E. Kemme [6] ran a series of 
experiments with a horizontal heat pipe (thus eliminating the gravity ef- 
fect) using liquid sodium and potassium in heat pipes which had wicks 
composed of rectangular channels. Kemme in some cases covered the channels 
with fine mesh screen and studied the effect of the screen on the per- 
formance. Ten heat pipes constructed of nickel were utilized. These tests 
demonstrated that the heat transfer characteristics of heat pipes with 
fixed outside dimensions can be changed significantly by variations in the 
wick structure. The practical temperature range over which such systems 
operate will be determined by the working fluid. Kemme [7] also investigated 


the entrainment limit and found that this phenomenon could limit heat transfer. 


16 


The entrainment limit occurs when the velocity in the vapor phase be- 
comes great enough to entrain droplets of liquid from the wick thus pre- 
venting the condensate from returning to the evaporator. Carnesale, et. 
al., [8] developed a theory of operation for heat pipes and experimentally 
investigated pipes which utilized water for the working fluid and packed 
monel beads for the wick. Carnesale also investigated the effect of the 
inclination of the pipe with respect to gravity in his experiments. 
Carnesale's group developed their own theoretical equations which are 
similar to equation (5). Their experimental results were in close agree- 
ment with their theoretical results. J. E. Deverall and E. W. Salmi [9) 
have tested a heat pipe in space and have found that the pipe operated 
normally. In addition, Deverall and Salmi ran an experiment by attaching 
a heat pipe to a vibration tester and operating the pipe as the vibration 
tester was in operation. They found that the vibrations had no adverse 
effect on operation and in some cases performance was improved because the 
vibration made the wick wet more rapidly and completely. D. M. Ernst [10] 
has performed further theoretical evaluation of Cotter's work [3] and has 
come up with some interesting modifications to Cotter's equation for maxi- 
mum heat flux (equation (5)). Ernst modified the boundary condition of the 
condenser and arrived at the following equation for the pressure drop in 
the vapor. 

AP (2) = 97/8 OR “1 (6) 
Samuel Katzoff [11] investigated the use of heat pipes and vapor chambers 
for thermal control in spacecraft. Katzoff discussed some of the practi- 
cal considerations concerning the use of heat pipes in space applications. 
Some of these considerations were wicking materials, working fluids, and 


the effect of contaminants such as non-condensable gases. The working 


Vy 


fluids which he investigated were alcohol, water, and a two-component 


fluid consisting of ethanol and methanol. 


Thesis Objectives 


The objectives of this thesis were to design and operate a heat 
pipe, and to determine what effect gravity and nucleate boiling in the 
wick have on heat pipe performance. 

The first objective was to determine the operational details of 
manufacturing a laboratory heat pipe for experimental use. Examples of 
the operational details sought were: (1) What difficulties were encounter- 
ed in the manufacture of the pipe? (2) Was the cleaning procedure correct? 
(3) Was the instrumentation adequate? (4) What factors should one keep 
in mind when designing a heat pipe? 

The effect of gravity appears only in the liquid pressure drop term 
in equation (5). The second objective was to determine if the effect of 
gravity on the liquid pressure drop is accurately represented and if 
gravity effects the operation of the heat pipe in any way not predicted 
by Cotter's theory. 

Cotter [3], Kemme [6], and others have stated that the formation of 
bubbles in the wick in the evaporator section of the pipe will limit the 
maximum heat flux that the pipe will transfer. In other words, the heat 
transfer mechanism of the pipe is evaporation and heat pipe operation 
fails when boiling occurs in the wick. Consequently, the final objective 
was to determine if boiling occurred during operation or at least at maxi- 
mum heat flux and to determine the effect of such boiling on the operation 


of the pipe. 
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II. EXPERIMENTAL EQUIPMENT AND PROCEDURES 
Equipment Design and Construction 

Once the application of a heat pipe has been decided upon, a suitable 
working fluid must be chosen. The working fluid must be capable of wet- 
ting the heat pipe material. It must also have high surface tension, high 
latent heat of vaporization and heat capacity, low viscosity, and a boil- 
ing point compatible with the operating temperature of the proposed heat 
pipe. In addition, the working fluid should, if possible, be readily 
available, easy to handle, and inexpensive. Water was selected for this 
application because it met most of the above criteria and the temperature 
range for this experiment was to be between 200°F and 300°F. Other 
liquids that have been used are sodium, potassium, and alcohol [6], [8], 
[il 

The heat pipe material should have high thermal conductivity and be 
compatible with the working fluid selected. Stainless steel, if properly 
cleaned, has good wetting characteristics with water. Stainless steel 
is difficult to machine because of its hardness. When stainless steel 
is used with water as a working fluid, great care must be taken to ensure 
that the steel is perfectly clean. 

The initial design called for a wick consisting of rectangular grooves 
or slots cut into the stainless steel pipe wall and covered with a fine 
mesh stainless steel screen. The effective capillary radius of the screen 
determines the maximum driving force for flow in the wick. The driving 
force is inversely proportional to the effective capillary radius, there- 
fore it is desirable to have the effective capillary radius as small as 
possible. For porous media, the pressure drop in the liquid is inversely 


proportional to the square of the capillary radius. It can be seen that in 


19 


trying to increase the driving force in the wick by decreasing the 
capillary radius it is possible to create a pressure drop in the 

liquid large enough to inhibit operation. This situation may be avoided 
by the use of covered channels. The driving force is still provided by 
the small capillary radius of the screen but the pressure drop is govern- 
ed by the size and shape of the channels. Since the pipe was to be re- 
quired to perform against gravity, a wire mesh having a very small ef- 
fective capillary radius was required to cover the channels. Very little 
information concerning the effective capillary radius of wire mesh was 
available. Katzoff [11] had run capillary rise tests on six wire mesh 
screens and from his data it appeared that a 300-400 mesh screen would 
lift a column of water about 16-19 inches. 

The final design for the heat pipe consisted of a 1.5 inch outside 
diameter stainless steel pipe having a wall thickness of 0.18 inches and 
a wicking material consisting of three layers of 400 mesh stainless steel 
screen. The wicking material was held in place against the inside of the 
pipe wall by a helical wire spring inserted inside the pipe. The large 
wall thickness was required to allow thermocouple wells to be drilled 
axially into the wall of the evaporator section. 

A series of tests were run to determine the effective capillary radius 
of the 400 mesh stainless steel screen selected and the height of a column 
of water that could be supported by the mesh. The screen selected was able 
to support a column of water 21 inches high. 

One of the theoretical assumptions was a large length to diameter ratio 
for the pipe. The length of the pipe had to be limited so that the pipe 
would operate in a vertical position with the evaporator uppermost, i.e., 
the length of the pipe could not be greater than the height of the column 


of water that the screen would support. In fact, it should be less. On 
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this basis, a length of twelve inches and a diameter of about one inch 
were selected. The pipe was cut 12.36 inches long to allow the 0.18 
inch thick end caps to fit in each end. The resulting inside length of 
the pipe, including the wick, was twelve inches. One end cap was fitted 
with a % inch stainless steel tube and vacuum valve to use for filling 
the pipe and was then welded in place. 

The reason for the use of the mesh wicking material rather than the 
rectangular grooves was because of the difficulty encountered in machin- 
ing the grooves in a stainless steel pipe. The grooves in the wick were 
to be cut 0.010 x 0.060 inches in a flat piece of stainless steel 0.125 
inches thick. The steel sheet was then to be formed into a cylindrical 
pipe by rolling and seam welding. There was no facility available local- 
ly that could roll the steel without damaging the grooves, so some con- 
sideration was given to designing a pipe of square cross-section. This 
scheme was abandoned because of difficulties in machining the small 
grooves in the stainless steel. Some of the difficulties encountered were 
a high rate of blade breakage and mechanical inaccuracies exemplified by 
slight vibrations of the milling machine which caused the blades to cut 
through or bend the metal between the grooves. The above-mentioned diffi- 
culties caused the construction of the pipe to fall behind schedule. 

The final design created a large pressure drop in the liquid phase 
and resulted in the maximum theoretical heat transfer being reduced by a 
factor of ten or from about 6800 BTU/hr to 680 BTU/Hr. In order to mea- 
sure the pressure in the liquid phase a series of six pressure taps were 
drilled into the pipe. Three taps were placed at each end of the pipe and 
spaced 0.7, 1.7, and 2.7 inches from the ends, respectively. 


The requirement to measure the temperature distribution in the pipe 
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Figure 2. Thermocouple Positions in Evaporator Wall 


22 


wall necessitated placing thermocouples in the wall of the evaporator. 
Six holes 0.040 inches in diameter were drilled an axial distance of 

one inch into the pipe wall. The location of these holes is shown in 
Figure 2. Much difficulty was encountered in drilling the holes for the 
thermocouples. Countless drills were broken and the drills that did not 
break had a tendency to bend. In one instance the drill penetrated the 
interior of the pipe and the hole had to be welded shut. In two other 
cases the drills came out of the outer surface of the pipe and these holes 
had to be closed also. 

After the pressure taps and wicking material had been installed in 
the pipe and prior to sealing the remaining open end, the pipe was first 
thoroughly degreased with acetone and dried by a heat lamp. After drying, 
the pipe was bright dipped ten to fifteen minutes in an aqueous solution 
of 12% (by volume) nitric acid at 160°F. The pipe was then rinsed three 
times in distilled water and the remaining end cap was welded into place 
and the holes for the thermocouples were drilled. The heat pipe was then 
heated to about 500°F and outgassed with a vacuum pump for twelve hours. 
At this point it was discovered that the vacuum pump in operation was 
faulty and oil was getting into the heat pipe. The cleaning process was 
then repeated with some modification. Since the pipe had been sealed and 
the only access to its interior was through the filling valve and the 
pressure taps which were covered on the inside by the wick, the acetone, 
acid, and distilled water each had to be forced through the pipe. The 
pipe was then heated and outgassed again using a different vacuum pump. 

A cooling coil of 3/8 inch copper tubing was then tightly wrapped and 
soldered to one end of the pipe to form a water jacket for cooling the 


condenser. The length of the condenser section was 4-3/4 inches. 
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The heating element for the pipe was THERMOCOAX wire and SAUEREISEN 
electrical cement. The heater wire consisted of an inner conductor of 
nickel-chromium alloy covered with an electrical insulating layer of 
compressed magnesium oxide. The layer of magnesium oxide was covered 
with a sheath of inconel. The wire was tightly wrapped around the evapora- 
tor section and then was covered with the electrical cement. The purpose 
of the electrical cement was to hold the heater wire in place and to pro- 
vide a more uniform heat distribution over the surface of the evaporator. 

The heat pipe was filled with water prior to the installation of the 
heating element. The pipe was filled by distilling 7.5 ml of water from 
a flask and condensing the vapor in the heat pipe. This amount of water 
was calculated to be sufficient to completely fill the wick and to have 
a slight excess of water remaining. With the heat pipe filled and ready 
to operate, the instrumentation was installed prior to operation. The 
thermocouples were installed in the pipe wall and along the side of the 
pipe. The thermocouples used were 40 mil copper-constantan thermocouples 
referenced against an ice junction. Manometers were connected to the 
pressure taps to measure the pressure drop in the wick. Copper-constantan 
thermocouples were also installed at the entrance and exit of the conden- 
ser coil to measure the temperature change of the cooling water. 

A seven inch length of 4 inch round stock was attached to the pipe so 
that the pipe could be mounted and rotated in a vertical plane about the 
horizontal axis of the arm. The pipe was then mounted and wrapped with 
asbestos until a layer of insulation about 3/4 inch thick covered the pipe 
completely except for the condenser coil. 

Power leads were then attached to a variac which was connected to a 
115 volt alternating current power supply. An a-c wattmeter and ammeter 


were installed between the variac and the heater element to measure the 
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power supplied and to monitor the current. The requirement to monitor 
the current was necessary because the heater element wire would only 


carry seventeen amperes without melting. 


Operating Procedures 


The first operation of the pipe was in the vertical position with 
the evaporator at the lower end. Power was supplied at the rate of fifty 
watts for about fifteen minutes until condensation drops began to form in 
the tubes leading to the manometers. Further operation caused the pres- 
sure in the pipe to build up and to force the water from the pipe into 
the tubes to the manometers resulting in the pipe drying out. Consequent- 
ly, at this point operation was stopped. The pipe was then emptied and re- 
filled with 25 ml of water. This amount of water was sufficient to fill 
the wick and to fill the vapor chamber to a height of one and one-half 
inches from the bottom with the pipe in a vertical position. Five of the 
six pressure taps were then sealed and the remaining tap on the condenser 
end of the pipe was connected to a mercury manometer so that the pressure 
inside the pipe could be measured. The pipe was then returned to the pre- 
vious vertical position with the evaporator in the lowest position and 
heat was applied in increments of fifty watts over a range from zero to 
seven hundred watts. At each increment the heat input was held constant 
until the temperatures in the pipe wall stabilized and a series of readings 
were taken of all temperatures and the internal pressure of the pipe. 
Noise inside the pipe was monitored with a stethoscope and nucleation 
could be heard at a heat input of fifty watts and above. The pipe was 
then rotated to the horizontal position and heat was again applied in 
various intervals from fifty to four hundred watts. At each interval the 


temperatures were allowed to stabilize and a series of readings were taken. 
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The temperature did not appear to stabilize at a power input of four 
hundred watts and this was the final data point taken. The value of 
maximum heat flux was determined by plotting the axial temperature distri- 
bution and observing the rate of heat input at which the temperatures do 
not stabilize. The maximum heat flux was taken as the maximum heat input 
at which thermal equilibrium was attained. The pipe was then rotated to 
a position 60 degrees from the vertical with the evaporator higher than 
the condenser and the power was increased in ten watt intervals beginning 
with twenty watts until the limit was reached. The pipe was then rotated 
to positions 45 and 0 degrees from the vertical with the evaporator upper- 
most and the process repeated in each position until the point of maximum 
heat flux had been reached. Noise level was also monitored at each posi- 
tion of the pipe, but no nucleation or boiling noises could be detected 
in the 60, 45, or O degree positions. Boiling could be heard in the 
horizontal position. 

The rate of flow of cooling water through the condenser was measured 
at each heat pipe position by means of a stop watch and a weighing tank. 

It was noticed that the heat pipe seemed to have a small leak. This 
was evidenced by the fact that the pipe would lose part of its vacuum 
if left overnight or longer. The pressure in the pipe at room tempera- 
ture should be the pressure of saturated water vapor at that given tempera- 
ture and should be lower than atmospheric pressure. When the pipe had 
cooled from a series of runs, the pressure would be the saturation pres- 
sure at the temperature of the pipe. The pressure would increase over- 
night, indicating that air had leaked into the pipe. To overcome this 
problem, or at least to get a common starting point for each series of 
runs, the pipe was operated in the vertical position, evaporator at the 


lower end, at a heat input of fifty watts for fifteen minutes with the 
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loading valve open. 

The reason for removing air from the pipe was because non-condensa- 
ble gases collect in the condenser during operation and reduce the volume 
of the condenser that the vapor can enter. Consequently, the wick area 
available for condensation is reduced because the vapor cannot reach all 


of the wick in the condenser due to the non-condensable gases present. 


Measurement Techniques 

In attempting to adequately determine the performance of an experi- 
mental heat pipe, the parameters to be monitored are the heat fluxes into 
and out of the pipe, the axial temperature distribution, and the pressure 
distribution. It is sometimes convenient to be able to watch or listen 
to the action inside the pipe during operation. 

The heat fluxes into and out of the pipe were measured in two dif- 
ferent manners. The heat into the pipe was determined by measuring the 
temperature at various points in the pipe wall and using the resulting 
temperature distribution in the wall and the equations for steady-state 
conduction assuming one-dimensional radial heat transfer to determine 
the heat flux. The heat out of the pipe was determined from the differ- 
ence in temperature and the mass rate of flow of the cooling water through 
the condenser. 

The axial temperature distribution was monitored by placing two 
thermocouples along the section of the heat pipe between the evaporator 
and the condenser. One thermocouple was placed on the outside wall of 
the pipe two inches from the evaporator section and the second thermo- 
couple was placed on the outside pipe wall four inches from the evapora- 
tor section and about % inch from the condenser. 

The pipe was not constructed to accommodate visual observation, how- 
ever it was possible to monitor the noise generated during operation by 


use of a stethoscope. 
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III. DISCUSSION OF RESULTS 

Theoretical Calculations 

The theoretical results were obtained by solving equation (5) for 
the heat flux. The properties of the fluid were taken to be the proper- 
ties of water at 212°F and one atmosphere of pressure. Pipe and wick 
dimensions were the actual dimensions of the pipe and wick. The contact 
angle, 0, was assumed to be close to zero so that cos @ is approximately 
equal to one. The wick void fraction was approximated by dividing the 
difference in weight of a sample of the wire mesh and the weight of a 
solid piece of stainless steel the same size by the weight of the solid 
piece of stainless steel. The screen permeability factor, b, was assumed 
to be twenty. 


Equation (5) is of quadratic form which may be readily solved for Q 


7 2 
Q = -B Ps B -4 AC G7.) 


2 A 


Where A is the coefficient of cap B is the coefficient of Q, and C is 


giving: 


composed of the remaining terms. As it turns out, very small values of 
the effective capillary radius, oe tend to make B so large that the 
equation yeilds either a very large negative value for Q or a trival solu- 
tion. Neither of these solutions is meaningful. To avoid this situation, 
equation (5) may be solved for maximum heat flow by assuming a value for 

Q and solving the equation by some iterative technique. Figure 3 shows 
the values of Q versus @. The theoretical values of Q were obtained 

by solving equation (5) by the Newton-Raphson method. A sample calcula- 


tion is given in Appendix I. 
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Experimental Results 

Most of the experimental results are shown in Figures 3 through 8. 
Figure 3 shows the theoretical and experimental values of Q versus angle 
of inclination, @. Figures 4 through 8 show the temperature distribution 
along the length of the pipe as a function of heat input at various 
angles of inclination, §@. 

The heat pipe did not operate as the theory predicted. For angles 
of inclination less than 90 degrees, the experimental results were much 
lower than the theoretical values. There are three possible reasons for 
this behavior. First, the pressure drop in the wick was much greater than 
anticipated or second, the working fluid was not wetting the wicking 
material so that capillary action was not taking place. Either of the 
two above phenomena would be detrimental to operation resulting in in- 
sufficient condensate flow to the evaporator for sustained operation 
except at much lower heat fluxes than those predicted. The remaining 
reason for the. large discrepancy between theoretical and experimental 
values for maximum heat flux could also be attributed to not knowing the 
exact positions of the thermocouples in the heat pipe wall. The error 
in Q that is introduced by small inaccuracies in the thermocouple posi- 
tions can be on the order of 100% (Appendix II). 

At an angle of inclination of 90 degrees and greater, the device was 
acting as a reflux boiler, causing the maximum heat flux to be much great- 
er than heat pipe theory would predict. 

Failure to operate close to the theoretical predictions can be at- 
tributed to manufacturing and cleaning errors which were not apparent until 


during or after operation. 
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Examination of Pipe After Operation 

After all the experimental runs had been completed and data collect- 
ed, the heat pipe was disassembled and one end of the pipe was cut off. 
The wick was removed and inspected. Inspection of the wick revealed that 
the wick and retaining spring had been partially decomposed and oxidized. 
The cause for this oxidation is not known, however it is entirely pos- 
sible that all of the acid was not removed from the pipe when it was 
cleaned the second time and this small amount of remaining acid could 
have caused the wick to decompose partially when the pipe was subsequently 
heated and outgassed. The presence of oxidation products in the wick 
would have had an adverse effect on the wetting of the wick and on flow 
through the wick. This non-wetting would promote nucleate boiling and 


boiling noise was heard during some of the runs. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 

Conclusions 

1. The results of the experimental work were inadequate to lead to 
any definite conclusions concerning the effect of gravity or nucleate 
boiling on the operation of a heat pipe. 

2. Great care has to be taken in the manufacture, cleaning, and 
filling of heat pipes if they are to operate properly. 

3. Position of thermocouples must be octane known in order to 
accurately determine the heat flux through the evaporator wall. 

4. The pressure distribution in the wick cannot be measured by any 
method that disturbs the integrity of the pipe. 

5. Prior to designing a heat pipe, one should obtain all the perti- 
nent information about the wick by experimental means if it is not pub- 


lished or if its validity is in doubt. 


Recommendations 

1. First, consider the measurement of heat flux into the evapora- 
tor. Measurement of power into the heating element does not accurately 
give the heat flux into the evaporator so some direct means of measure- 
ment is required. One method would be to machine small grooves in the 
outer and inner walls of the evaporator, place thermocouples in the 
grooves, accurately determine their positions, cover the thermocouples and 
fill in the grooves with the heat pipe material, and polish the surface 
so that it is even with the adjacent material. Another method would be 
to measure the heat flow directly. Geonetics Corporation of San Diego 
manufactures a device called a Thermal Flux Transducer that could possibly 
be placed on the evaporator surface between the pipe and the heating ele- 


ment and would measure the heat flux directly. 
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2. Measurement of axial temperatures in the vapor chamber can be 
accomplished by inserting a small tube containing the desired number 
of thermocouples. This could be done with a minimum of difficulty using 
a CONAXmultiple seal gland. 

3. Small pressure transducers installed in the wali of the pipe 
should serve to measure the pressure distribution in the wick while 
maintaining the Titec xivesy of the interior of the pipe. 

4, The following method is recommended for determining the proper- 
ties of a porous wicking material prior to installation in a heat pipe. 
Three important properties of the wick material are the capillary radius 
ro the permeability factor b, and the screen void factor, € . The 
value of r. for a particular wick may be determined by a simple capil- 
lary rise experiment. The factors b and © appear only in the term for 
the pressure drop in the wick and may be grouped as b/€. A value for 
this term may be determined as follows: 

Take a section of wick material of known constant length and 
cross-sectional area, pass a constant, steady flow of the working fluid 
through the material and measure the pressure drop. The Liquid pres- 


sure drop term in equation (5) is 


Y_2 b 
2 (Bn -Bo) r = 


rtj{o 


Now Q/L = m= mass flow rate, 
oe 2 ‘ ‘ 
2 widiwis' -~R') represents the cross-sectional area, Z is the length of 


the sample, \) is known and r, can be found. Therefore 


a x area 
= ae 8 
b/e = AP, LIQ \) ® length St 
5. The cleaning procedure should be thoroughly investigated and 


tested prior to using. 
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6. If at all possible, the heat pipe should be designed so that 
visual observations can be made of the interior of the pipe during 
operation. 

hes The condenser end of the pipe should be flanged so that the 
pipe can be dismantled for examination after operation or for recleaning 
if necessary. 

8. Build a heat pipe of rectangular cross-section with a wick on 
either one or three sides. Heat would be added or removed from the sides 
covered by the wick. The fourth side would be enclosed by either two 
or three parallel layers of glass, each sheet of glass separated from 
the other by a small distance and this space between the sheets of glass 
should be evacuated. Prior to operation, the entire pipe should be 
uniformly heated at the operating temperature to insure that there will 
be no condensation on the glass. This device could also be instrumented 


for pressure and temperature measurement. 
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APPENDIX I 
Sample Theoretical Calculations 
The solution to Cotter's pressure drop equation (5) for maximum 
theoretical heat flux was determined by the Newton-Raphson method for 
solving a quadratic equation [12]. The Newton-Raphson method is an 


iterative technique and in its basic form is given by 


X41 x, -£ (x, )/£ (x, ) 


or in terms of Q 


Q ext 7 24-F(Q,)/£'(Q,) 


next 

Let A = (1-4/7 *)/@ RtL?) 
B= (bY 2)/(2 TH(R,*-R*) € x_“L) 
C= mm g Z cos O/g-2 ¥ cos O/r 


Cotter's equation may now be written 
AQ?+BQ4C=0 
Let 
2 
£(Q) = AQ +BQtC 
and 
£'(Q) = 2 AQ4B 


The following calculations are for water and steam at 212°F. 


R = 0.570/12 ft 

WwW 

R = 0.561/12 ft 
assume 

cos 9g = l 

cos 9 = l 

b = 20 

therefore 


A = .01395 Lb=sec”/BTU =£t" 
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2 


B = 1032.2 lb-sec/BTU-ft 


C -52.718 lb/ft? 

There are two roots, one positive and one negative, to the quadratic 
equation. The negative root is meaningless, therefore the positive root 
is the desired solution to the equation. The positive root will be a 
small positive number, therefore assume a value 


Q, = 0.12 


as a starting point. 


Q * .12-71.12 Ze .06 
pest 1032.2 
Since ~12-.06 O00] 
| see | ? 
set QQ = Q erxee = .06 


and iterate again following the above procedure until 


fee |e 
gy 





The final value of Qa? the desired root is 


Q, “= .051 BTU/sec 
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APPENDIX II 
Error Analysis 
One major source of error in determining the maximum capacity of 
a heat pipe lies in the inaccuracies in determining the heat flux at 
which operation breaks down. Direct measurement of power into the device 
by means of a wattmeter or other instrument is inaccurate due to inherent 
heat losses involved. An attempt was made to measure the heat flux by 
measuring the temperature at various points in the pipe wall and comput- 
ing the heat flux by means of the equation for steady state conduction 
in one dimension. The following error analysis is based on the heat 
transfer equation for one dimensional, steady state heat transfer in a 
cylindrical system which is just 
Q=2 Wk2(2,-T,)/In(r(/r,) 
for heat flow into the pipe. The first four terms on the right side of 
the equation may be grouped into one constant term so that 
Q = K(T,-T, )/in(r,/r, ) 


Now setting QT =T -T, 


and taking the derivative 


3Q= K SAT) + i SK . K AAT S(r6/r,) 2 
In(r r,) In(r, r,) er fiw) 


So-Sar , Sx - 68!" 
Q AT K (=./7;) In/ ral, ) 
Now K is a constant term and may be assumed to be known accurately with 


respect to the other terms so that 


41 


Also, since the desired value for $ Q/Q is its maximum value, the 


negative term is replaced by its absolute value; therefore: 


=$o7 + Cn) 


A (r,/*,)10 (r/r.) 


4 
Foo) Feo) 


Now assume the following values: 


Q = 50 watts = 170.65 BTU/hr 


ro = 0.720 inches 
vr, = 0.625 inches 
K = 2Wkz = 15.71 BTU/hr-°F 


Therefore: 


(170.65) 1n(. 720/625) 
at = 15.71 


= 1,59°F 

Now the error in measuring the temperature was due to two causes. The 
first was the inaccuracies in the digital voltmeter and they were ap- 
proximately +2 Av. The second cause of error arises from the method of 
calibrating the thermcouples. The thermocouples were calibrated in a 
steam bath at 212°F and the error was assumed to be constant throughout 
the temperature range measured. The error in measurement due to the 
above assumption was estimated to be approximately +8 AAv. The total 
temperature error was 8 +2 or 10 AAv. An error of 10 JA v is equi- 


valent to a temperature error of approximately 0.4°F for copper-constantan 


thermocouples. This gives an estimate of 


54) = Oa) = .25 


AT 1.59 
The remaining error was one of thermocouple position error and was based 
on two premises. The first premise was that the holes for the thermo- 
couples were not true and the second was that the location of the thermo- 


couple junction in the hole was not known. 
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For purposes of error analysis the assumption was made that the 
end of the drill bit could deflect one bit diameter between the en- 
trance and the end of the hole, that is, the center of the entrance 
hole and the center of the bottom of the thermocouple well either 
are not the same radial distance from the center of the pipe, or have 
some angular displacement between them, or both. The average size of 
the thermocouple bead was 0.018 inches in diameter, including insulation. 
It was assumed that the copper-constantan junction was in the center of 
the bead. Under the above conditions it would be possible for a given 
thermocouple to have a maximum error in radial position, Ar, of 0.020+ 
0.040-0.018/2=0.051 inches. An error of this magnitude would cause a 
maximum error in the ratios of two different radii, S (x /r,); of 0.169 
when ee and r. are the values given above. 


This gives an estimate of 


S eeu = 169 = 1.04 
(©,/*;)48 (ro/r,) (1 US Aled 5) 


Therefore for maximum error $ @/Q is: 


Sao FS i207. 
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